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ABSTRACT 

Hydrogen-rich Type Il-PIateau supernovae exhibit correlations between the plateau luminosity Lpi, the nickel 
mass Mni, the explosion energy iJgxp, and the ejecta mass Mej. Using our global, self-consistent, multi-band 
model of nearby well-observed supernovae, we find that the covariances of these quantities are strong and that 
the confidence ellipsoids are oriented in the direction of the correlations, which reduces their significance. By 
proper treatment of the covariance matrix of the model, we discover a significant intrinsic width to the correla¬ 
tions between Lpi, ii^exp and Mni, where the uncertainties due to the distance and the extinction dominate. Eor 
fixed i?exp, the spread in Mni is about 0.25 dex, which we attribute to the differences in the progenitor internal 
structure. We argue that the effects of incomplete 7-ray trapping are not important in our sample. Similarly, 
the physics of the Type Il-Plateau supernova light curves leads to inherently degenerate estimates of E^exp and 
Mej, which makes their observed correlation weak. Ignoring the covariances of supernova parameters or the 
intrinsic width of the correlations causes significant biases in the slopes of the fitted relations. Our results 
imply that Type Il-Plateau supernova explosions are not described by a single physical parameter or a simple 
one-dimensional trajectory through the parameter space, but instead reflect the diversity of the core and surface 
properties of their progenitors. We discuss the implications for the physics of the explosion mechanism and 
possible future observational constraints. 

Keywords: Methods: statistical — stars: distances — supernovae: general 


1. INTRODUCTION 

The observed light curves and expansion velocities of 
hydrogen-rich Type Il-Plateau supernovae can be used to 
infer the properties of the explosions and progenitor stars. 
Specifically, the duration and luminosity of the optically-thick 
“plateau” phase of nearly constant bolometric luminosity is 
primarily set by the explosion energy Ep-^r. and ejecta mass 
Mej (e.g. lArn^ll98(]HKasen & Woosley^^2009^ . The subse¬ 
quent nearly exponential fading is powered by the thermal- 
ization of radioactiye fission products of ®®Ni and the lu- 
minosity is th us proportional to the nickel mass Mni (e.g. 
iHamuyl120031) . Patterns in the distributions and the correla¬ 
tions between these quantities can guide the stellar eyolution 
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ample, it has been proposed that the supernoya explosion 
energy is proportional to t he ejecta mass and therefore also 
the p r ogenitor mass (e .g. iHamuyl 120031 : lUtrobin & Chugai 
l2009t iPoznansO l2013h . Howeyer, low-energy explosions 
might be an exception signaling significant fallback in mas- 
progenitors (e.g. IZampieri et akl l2003t iPastorello et al.l 
At the same time, direct progen- 
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itor detections show significant scatte r but little co rrelation 
between the progenitor mass and Mni (ISmartlll2009h . 

Naturally, not all supernoya parameters can be inferred in¬ 
dependently. Eor example, quantities based on bolometric 
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luminosity such as Mni or the plateau luminosity Lpi are 
plagued by systematic uncertainties in the distance or ex¬ 
tinction. Similarly, a simultaneous change in seye ral pa- 


rameters can result in nearly ident 
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ical light curyes ( 
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Here, we eyaluate whether or not the systematic uncer¬ 
tainties and parameter coyariances influence the significance 
of the correlations between parameter estimates for Type II- 
Plateau supernoyae. We focus on Lpi, Mni, Lexp, and Mej. 
To this end, we employ the self-consistent global model of 
nearby we ll-obseryed Type Il-Plat eau supernoyae that we de- 
yeloped in iPeicha & Prietol (120151 hereafter PP15), which si¬ 
multaneously fits multi-band light curyes and expansion ye- 
locities and proyides consistent distances, reddenings, bolo¬ 
metric luminosities and, most importantly for the present pur¬ 
poses, their coyariances. In Section |2] we describe how we 
estimate the supernoya parameters, the database of supernoya 
obseryations and the estimates of uncertainties. In Section [3 
we inyestigate the systematic uncertainty due to distance and 
discoyer a significant intrinsic width of the Lexp-MNi corre¬ 
lation. In Section m we address the significance of the Mgj- 
Lexp correlation. In Section |5] we discuss the astrophysical 
implications of our findings. 

2. SUPERNOVA PARAMETERS DERIVED FROM OBSERVATIONS 

We calculate the bolometric luminosity of the optically- 
thick plateau phase Lpi as 

Upi — Lboi(fo “f Afpi), (1) 

where Lboi is the bolometric light curye obtained by integrat¬ 
ing the spectral energy distribution from about 0.19 to 2.2 /rm 
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and extrapolating the Rayleigh-Jeans tail dPPlSh . and to is 
the zero point of our model fits, which coincides with the 
explosion epoch for the purposes of this paper. The lumi¬ 
nosity is evaluated at a f ixed interval A fni after the explo¬ 
sion epoch to- Following iHamuvl (120031) . our fiducial choice 
is Afpi = 50 days, but we will investigate the sensitivity of 
some of our results to Afpi. We estimate the nickel m ass Mn; 
from t he exponential decay tail of the light curve after lHamuvl 
(l2003l) as 

ATni = 7.866x10 ^^Lboi(fo+^iNi) exp ^^ Mq, 

( 2 ) 

where r = 111.26days. Our fiducial choice for the time 
elapsed after explosion, where we estimate Mni, is A^Ni = 
200days. The units on Lboi are ergss“^. 

We estimate the explosion energy E'exp and ejected mass 
Mej using linear relations of the form 

where b = {My, logfp, logu), and My is the absolute mag¬ 
nitude in the V band, tp is the duration of the optically- 
thick plateau phase measured at the midpoint of the drop 
to the exponential decay phas^E and v is the expansion ve¬ 
locity of the photosphere in the units of 1000km s“^ com¬ 
monly measured on the Fe II 5169 A line. My and v are 
evaluated at the midpoint of the plateau, corresponding to 
time to + tp/2. The coefficient vectors a and /3 are typ¬ 
ically obtained either from analytic models of supernova 
light c urves (e.g. lArnetlll980l : rPopovll9^lKasen & WooslevI 
120091) or from the fits t o the simulated light curves and ex¬ 
pansion velocities (e.g. lOtvinova & NadezhuH I1983L 119851 
iKasen & WooslevI 1200^. Hp e. we use coefficients from 
the analytic model of IPopovI (119931) . a. = (0.4,4.0, 5.0), 

(3 = (0.4,4.0,3.0), = -3.311, and 77. 1 = - 2.089, 

and the simulations of lLitvinova & NadezhinI (119851) . a = 
(0.135,2.34.3.13), 0 = (0.234,2.91,1.96) = -3.205, 

and 77ej = —1.829. iLitvinova & NadezhinI (1 19^ claim that 
Equations (I3-0|i can reproduce their numerical results to about 
30% and t heir results have been commonly used in the liter¬ 
ature (e.g.lElmhamdi et ^120031 [Hamuy|[2003l [Hendry et alJ 
I2006t iBose et al.ll2013l) .~We emphasize that the point of this 
exercise is not to compete with detailed modeling employing 
sophisticated codes, but to illustrate the limitations posed by 
the physics of the supernoya light curyes to the estim ates of 
.Egxn and Mej - the analytic results for a and (3 of iPopovI 
(119931) already imply that the estimates of iSexp and Mej are 
highly correlated. We will show in Section]?] that similar de¬ 
generacy is present also in the more sophisticated models. 

To estimate the aboye para meters , we use the model and 
database of obseryations from lPPlSl The light curyes are de¬ 
scribed by a set of phenomenological parameters a obtained 
by least-squares fitting of the data, which also proyides the 
full coyariance matrix C® of our model. This allows us to 
properly propagate uncertainties in a to Lp\, Mni, i?exp, and 
Mej. In this work, we employ a model fit with obserya- 
tional uncertainties rescaled so that the final reduced is 

^ The plateau duration tp should not be confused with the time, when the 
plateau luminosity is measured, Atpp 


unity, which increases the yalues in C® relatiye to the un¬ 
adjusted fiB We use a subset of the obseryational sample 
of lPP15l that does not include SN2007od, SN2006bp, and 
SN2002hh due to reasons mentioned in 1PP15L In addition , 
we include a recent Type II-P SN2013am dZhang et al.l]2M4l) . 
which we fit using the publicly ayailable yersion of our fit¬ 
ting tool0. Eor SN2013am, we find an explosion epoch to = 
2456373.0±2.4, total reddening F;(B-1/) = 0.81±0.02, and 
a distance modulus 29.2 ± 0.3 mag. Eor eyery choice of Afpi 
and AfNi we select a subset of supernoyae with data before 
and after these dates to preyent extrapolation of the model. 
Our fiducial choice of Afpi and Atya thus leayes us with 19 
supernoyae, which satisfy these constraints. 

We calculate the coyariance matrix of f = 

log(Lpi, Mni, E’exp, ?7ej) using the standard procedure for 
uncertainty propagation 



The confidence ellipsoid for / is a quadratic equation in the 
offsets 5f 

Ax^ = (6) 

where Ax^ depends on the desired confidence leyel and the 
number of yariables. In the subsequent discussion, we will 
exclusiyely focus on pairs drawn from f, (/j, fj), their 2x2 
coyariance submatrix and the 6 8.3% confidence leyel, 

where Ax^ ~ 2.30 (iPress et al.lll992l p. 697). Eor two pa¬ 
rameters at a time. Equation (]6]l can be solyed to obtain the 
positions on the confidence ellipsoid by transforming the pair 
of yariables (5f = {Sfi, Sfj) to polar coordinates, yarying the 
polar angle, and solying for the radial distance at each po¬ 
lar angle. Note that the off-diagonal elements of can be 
non-zero eyen if is diagonal. In our coyariance matrix, 
we properly take into account the uncertainties in to and their 
associated coyariances with the distance and plateau duration. 

We model the dependencies between indiyidual compo¬ 
nents of f = log(Lpi, Mni, ii^exp, Afej) with a straight 
line allowing for an intrinsic width S. We use the like¬ 
lihood function of iHogg et al.l (1201 Oh . which assumes that 
the obseryations are offset from a linear relation by a Gaus¬ 
sian described by the coyariance matrix conyolyed 

with a Gaussian intrinsic scatter with standard deyiation E, 
which is perpendicular to the linear relationship. We ob- 
tain the confidence interyals w ith the MCMC sampler emcee 
dEoreman-Mackey et al.|]20T3]) . We will quote median of the 
distribution for our best-fit parameters, and 16 and 84 per¬ 
centiles as their confidence interyals. 

We quantify the significance of the correlations in the data 
using two approaches. Eirst, we calculate the Bayes factor B 
of the linear fit relatiye to a model, which assumes no correla¬ 
tion between the two yariables. Specifically, we eyaluate 

n _ 71- J-71-/2 V > _ 

“ max[P(0 = O),P(0 = 7r/2)]’ 

^ We repeated the analysis with unadjusted covariance matrix and found 
that line slopes, intercepts and their uncertainties remain unchanged. The 
intrinsic widths are about 20% higher when using the unadjusted full covari¬ 
ance matrix. The Bayes factors B change by only about 20%, but K increases 
in all cases by a factor of ~ 2.5. Finally, the unadjusted covariance matrix 
will yield smaller confidence ellipsoids in the Figures. 

^ http: / /v!V!vi .astro . prince ton. edu/ ~pe j cha/iip/ 
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where 6 is the angle between the line and the x axis and P{0) 
is the likelihood marginalized over the line intercept and in¬ 
trinsic width together with their priors; we assume flat prior 
in 6. Taking the maximum of = 0) and P{9 = tt/2) 
in the denominator of Equation O ensures that B ni 1 for 
uncorrelated data when the scatter in the x and y directions 
differs. Equation (|7]i is also insensitive to rescaling the uncer¬ 
tainties of the data, because these enter in the same way both 
in the numerator and the denominator. According to I Jeffrey^ 
(119831) . B > 10^/^ implies that the support for the fitted line 
is “substantial”, and “decisive” if i? > 10 ^. 

The employed likelihood model is only concerned with dis¬ 
placements of the observations perpendicular to the fitted line 
and provides no in formation on the distribution along the line 
dHogg et alj|2010l) . To give a quantitative measure of the dy¬ 
namic range of the data along the best-fit line, we project the 
data and their confidence ellipsoids on the best-fit line and cal¬ 
culate the weighted standard deviation of the data along the 
best-fit line V and the median uncertainty along the best-fit 
line W. We define 

In the case of one-dimensional data, ii is a measure of the 
intrinsic scatter, for example, R = 3 would imply approxi¬ 
mately 3(7 significance of the intrinsic scatter. Here, the in¬ 
formation provided by R is complementary to the fitted slope 
and its uncertainty. Very large values of R together with small 
relative uncertainty on the slope and high B imply a strong 
correlation. If the confidence ellipsoids are oriented along the 
observed correlation and the off-diagonal elements of are 
ignored, R will be artificially higher and the observed correla¬ 
tion will appear stronger. Unlike B, the parameter R depends 
on the absolute values of data uncertainties. 

3. INTRINSIC SCATTER IN THE Lpi, Eexp AND Mni 
CORRELATION 

Significant off-diagonal terms in occur when the un¬ 
certainties are dominated by a single systematic uncertainty 
Csyst, typically in the distance or the extinctioifl Quanti¬ 
ties linearly proportional to Lboi such as Lpi and Mni are 
particularly susceptible. In other words, a bias in the dis¬ 
tance will move Lpi and Mni in the same direction simul¬ 
taneously by the same amount, introducing a covariance in 
these two parameters. Schematically, the covariance matrix 
of log(Lpi,MNi) is 


QlogLpi.logMNi 


^pl 3” ^syst ^syst \ 
'^syst ^tail '^syst / 


(9) 


where Upi and Utaii the uncertainties in the observed mag¬ 
nitudes during the plateau and the exponential decay tail, re¬ 
spectively. Usually, CTpijCTtaii crgyst and the confidence el¬ 
lipsoid is strongly elongated in the direction of the systematic 
uncertainty. If log Lpi and log Mni are perfectly correlated, 
neglecting the off-diagonal terms in Equation (|9]l will imply 
a value of R that is twice the true value. Even if the full co- 
variance matrix is not available, an approximate covariance 
matrix similar to Equation (0 can be constructed to properly 
visualize the confidence ellipsoids. 


^ Our distance estimates are compared t o previ ous results and other tech¬ 
niques (Cepheids, Type la supemovae) in IPPl.'il and yield generally good 
agreement. Similarly good agreement is obtained for reddenings. 


In the left panel of Eigure [T] we show the estimates of Lpi 
and Mni with uncertainties simply represented by the diag- 
onal terms of C^, as is commonly don e (e.g. iHamuvl 120031 
lAnderson et al.ll20T4t ISniro et al.ll2014l) . We would infer that 
there is a linear correlation between log Lpi and log Mni with 
a slope of 1.511 q} 7 and R = 4.2, implying a strong correla¬ 
tion. More importantly, considering only the diagonal terms 
gives a false impression that all the points are compatible with 
the best-fit line given their uncertainties. Quantitatively, there 
is no evidence for intrinsic scatter, with E = O.OSlQgg. 

When the confidence ellipsoids are properly included as in 
the right panel of Eigure [T] the picture changes. The confi¬ 
dence ellipsoids are significantly elongated, because of their 
mutual dependence on distances, as indicated by the arrovv 0 . 
The correlation is noticeably less significant, R = 3.0, al¬ 
though there is no doubt this correlation exists given the 
large dynamic range of the parameters. The Bayes factor is 
B 9 X 10^ implying strong support for the correlation. 
More importantly, we discover a statistically significant in¬ 
trinsic width of the relation E = 0 . 12 ^q q 2 > which implies 
a scatter of 0.2 dex in Mni for a fixed Lpi. Eurthermore, ne¬ 
glecting the off-diagonal terms or the i ntrinsic width of the re - 
lation can bias the inferred slope (e.g. iTremaine et al.ll200^ . 
Neglecting the off-diagonal terms increases the slope by about 
0.19 with a corresponding change in the intercept. Not ac¬ 
counting for the intrinsic scatter leads to slopes of l.Sd^g gg 
and for the full and diagonal covariance matrix, re¬ 

spectively. 

Since Lpi does not have an immediate physical interpreta¬ 
tion, we show Mni as a function of i n Eigu re |2] for the 
scaling relations of iLitvinova & NadezhinI (119851) and iPopo'^ 
(119931) The relative position of the majority of the data points 
remains unchanged when compared to the right panel of Eig¬ 
ure [T] which indicates that Lpi is a good proxy for E^xp- 
There are small differences between the two scaling rela¬ 
tions, but the relative p ositions of the majority of the points 
are unchanged. For the ILitvinova & NadezhinI (119851) coeffi¬ 
cients, we find that the Lexp-AfNi correlation is less signif¬ 
icant than Lpi-MNi with R = 2.8 and 3.7 for the full and 
diagonal covariance matrix, respectively. The Bayes factor is 
B « 1.7 X 10^ indicating strong correlation, but weaker than 
Lpi-MNi. The inferred intrinsic width orthogonal to the line 
is slightly higher than in the Lpi-MNi correlation but again 
statistically significant, E = O.ldlggg or 0.25 dex in M ni for 
fixed Lexp- The intrinsic width from the iPopo^ (119931) cali¬ 
bration is E = O.lOlgg®. 

The intrinsic width of the Lexp-AfNi correlation could 
be due to the 7 -ray trapping efficiencies Aj varying 
among supernovae with the same Lexp- Since the 
exponential decay luminosity is proportional to [1 — 
exp(—exp( —AfNi/r) (e.g. IChatzopoulos et al.) 
l2012tlNagv et al.ll2014l) . where AfNi is the time elapsed since 
the explosion, supernovae with significant 7 -ray leakage not 
only appea r fainter at any point of this phase but also de¬ 
cay faster (lAnderson et al.ll2014h . and their light curves di¬ 
verge from those of supernovae with full 7 -ray trapping over 
time. Since our sample contains supernovae with decay rates 
compatible with full trapping as evidenced by the exponential 

® Note that the unce rtainties in absolute magnitude and expansion velocity 
(e.g. IPoznansfl 1201311 should not be very conelated, unless the velocities 
were used for an estimate of the distance modulus, in which case there should 
be a signihcant correlation. 
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Figure 1. The correlation between the plateau luminosity Lpi (Eq. [T]) and the nickel mass Mni (Eq. (^). The best-fit relation is shown with solid line and the 
intrinsic width S is indicated with dashed lines. Left. The relation as usually presented, treating the pai'ameter uncertainties as if they are uncoiTelated. We show 
one-dimensional uncertainty projections and the correlation does not exhibit internal spread. Right: The full confidence ellipsoids for each supernova imply a 
significant internal spread in the correlation. The an'ow shows the direction of the covariance created by the distance uncertainty. 


Litvinova & Nadezhin (1985) calibration Popov (1993) calibration 




Figure 2. Nickel mass Mni as a function of explosion energy F^exp, with t he confidence ellipses properl y visualized. T he colors of t he individual supemovae 
are the same as in the right panel of FigurefJ We use the scaling relations of ILitvinova & NadezhinI ri985l left panel) and lPopo^j 19931 right panel). 




A?Ni [days] A?p, [days] 

Figure 3. Weighted standard deviation of the bolometric magnitude difference between the plateau and the exponential tail are shown with filled circles as 
function of (left panel) or At^i (right ^nel). There is little dependence on indicating that the different 7 -ray trapping efficiencies are not responsible 
for the intrinsic width seen in Figures^andlz] In other words, the exponential decay is nearly parallel for stars in our sample. Open triangles indicate the mean 
uncertainty of the individual measurements used. Only supernovae with data spanning Afpi and AfNi are used. 
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decay slopes dPPlSh . the 0.25 dex difference in the inferred 
Mni at A^Ni = 200 days should increase to about 0.7 dex at 
AfNi = 400 d, if the scatter is due to 7 -ray leakage in some 
objects. 

To test whether the late light curves of supernovae in our 
sample diverge with time due to incomplete 7 -ray trapping 
in some objects, we show the weighted standard deviation of 
the bolometric magnitude difference between the plateau and 
the exponential tail as a function of the time elapsed since 
explosion, AfNi, in the left panel of Figure[3] We see that the 
bolometric magnitude dispersion increases from 0.37 mag at 
AfNi = 200d to 0.45 mag at AfNi = 400 days, much less 
than what we would expect if some supernovae showed full 
trapping and some only partial. This means that the slopes of 
the exponential decay are very similar among our objects and 
are compatible with full 7 -ray trapping. 

For the sake of completeness, we test what is the impor¬ 
tance of when is the plateau luminosity determined. We show 
the weighted standard deviation of the bolometric magnitude 
difference between the plateau and the exponential tail but 
now as a function of Afpi in the right panel of Figure [3 
For small Afpi, the dispersion is relatively high, presumably 
due to differences in the properties of the shock-heated ejecta 
shortly after shock breakout, but for Afpi > 40 d the disper¬ 
sion remains approximately constant. We conclude that the 
intrinsic width of the Eexp-Myn correlation is robust with re¬ 
spect to when exactly the plateau and exponential decay tail 
luminosities are measured, and that it is not due to variations 
in the 7 -ray leakag^H 

4. THE CORRELATION BETWEEN M^j AND Sexp 

Correlated uncertainties in f = log(Lpi, Mni, it^exp, Afej) 
can also occur when some of the vectors dfi/dh are nearly 
parallel. Equations (|3® imply that the covariance of log M^j 
and log Flexp is 


a 


log Mej ,log TJexp 


\/a 


log Mej Jog Mej Clog Bexp Jog Bexp 


a ■ (3 

RM’ 


( 10 ) 


if the uncertainties in the individual components of b are ap¬ 
proximately the same. This is reasonable, because 5 to 10% 
uncertainties in the distance modulus, fp and v are expected. 
We obt ain a high correlation of 0.94 and 0.97 for the a and 
f3 from iLitvincwa & NadezhinI (Il985h and lPopo^ (Il993h . re¬ 
spectively. This degeneracy comes the physics of the Type II- 
Plateau supernova light curves, where higher kinetic energy 
makes the material transparent earlier, which can be com¬ 
pensated for by increasing ejecta mass to produce a pproxi¬ 
mately the same plateau duration and luminosity (e.g. Arnetj 
19801: iLitvinova & Nadezhinlll985l iKasen & Wooslev 120091 
Dessart et al.ll2010ll2013h . T^e degeneracies between param- 
et ers describing supe rnova light curves were also investigated 
bv lNagy et al.lOOl ^ using semi-analytic models. As a result, 
inferences of -Eexp and Mej will be highly correlated in any 
technique based on light curves and expansion velocities. 


^ An additional piece of anecdotal evidence against significant 7 -ray leak¬ 
age comes from comparing SN2013am and SN2005cs, which have neai'ly 
identical luminosities for the first 70 days. However, SN2013am has 
a noticeably shorter tp and a higher infen'ed Mni than SN2005cs. This 
implies that SN2013am has slightly smaller Eexp and significantly smaller 
Mej than SN2005cs (Fig.|4). If 7 -ray escape were important, we would ex¬ 
pect smaller infen'ed Mni and faster exponential decay in SN2013am than 
in SN2 005cs. Yet, the ex ponential decay slope is almost the same in both 
objects jZhang et alJ2Q14l) and SN2013am has higher inferred M^[. 


In Figure |4] we show Mej and E'exp as inferred for our sam¬ 
ple and the two scaling relations. The uncertainty ellipsoids 
are elongated along the correlation, as expected. We find 
significant system atic offsets b etween the scaling relation¬ 
ships, in particular lPopovI (119931) calibration produces smaller 
and more realistic Mej, and smaller F^exp- The relative po¬ 
sitions of the indi vidual supernovae remains u nchanged in 
most cases. For the lLitvinova & NadezhinI (Il985h calibration, 
we find a relatively high uncertainty for the slope 2.09'^q'^q 
and a statistically significant intrinsic width to the correlation, 
E = O.ld^p pg. We find R = 2.0 and the Bayes factor of 
B 200, which implies that this correlation is much weaker 
than Fip,yn-MNi, albeit the evid ence is still “decisive” in the 
cla ssification of Jeffrevsl (Il983h . The results are similar for 
the lPopo\^ (11993 ) calibration. 

There are two outliers to the Mej-F^exp correlation: 
SN1995ad and SN1980K. The Ty pe II-Linear SN1980K has 
well-constrained distance dPPlSl) and relatively good light 
curve covering the transition to the exponential decay. Re¬ 
cent investigations of large samples of hydrogen-rich super¬ 
novae indicate that there is a continuum of li ght curve shapes 
betwe en Type Il-Plateau and Type II-Line ar (lAnderson et al. 
120141 and lSander s et al.ll20r5l but see also lArcavi et al.ll2012 
and iFaran et al.ll2014t) . suggesting that a larger unbiased sam¬ 
ple of supernovae might fill the space between SN1980K and 
the rest of our sample. 

SN1995ad has distinctively shorter tp ~ 85 days than 
most s upernovae in our sample, as found also bv llnserra et al.l 
(120131) . which is res ponsible for t he out lying results. The 
explosion epoch of llnserra et al.l (1201 3h is about 12 days 
later than ours and their adopted distance modulus is 
about 0.3 mag closer (although the distances are compat¬ 
ible within their la limits) than what we obtained in 
|PP3 From radiation-hydrodynamical modeling llnserra et al.l 
(120131) found log(£’exp, TFej) ~ (0.3, 0.7) in the units of 
Equations (OJ-Bl), w hile we find (0.4, 0.3) and (0.9, 0.8) us¬ 
ing the relations of lPopovI (119931) and ILitvinova & NadezhinI 
(119851) . respectively. Without the outlying SN1995ad and 
SN1980K, the Mej-Flexp correlation exhibits a steeper slope 
of 2.601 q 3() (compatible with the result of lPoznanskill20 1 3l) . 
significantly smaller intrinsic width E = 0 . 051 qq 3 , and 
R = 1.9. The Bayes factor increases to i? Ri 1.8 x 10^, 
as expected when outliers are removed. 

Since we are not interested in the specific values of E^xp 
and Mej, but rather in the process of their estimation, the 
results of the analytic scaling relations of Equations (13-01l 
are only approximate and must be confronted with more 
detailed models. In Eigure |5] we show tp and v of the 
radiation-hydrodynamic models of explosions o f non-ro tating 
red supergiant progenitors of iDessart et al.l (120101) with 
information on E^xp and Mej encoded in symbols and their 
colors. The models do not include heating by radioactive 
nickel and thus underestimate tp. The coverage of the 
parameter space by the theoretical models is not uniform, 
and in particular some observed supernovae fall in the areas 
where there are no models, making interpolation in these 
theoretical models to estimate the supernova parameters 
difficult. More importantly, the bulk of the observed su¬ 
pernovae cluster at tp ^ 120 days and v « 3500km s“^, 
and these parameters can be explained by a range of models 
going from low explosion energy (F^exp = 0.1 x 10 ®^ ergs) 
and low ejecta mass (Mej ss SMq) to normal explosion 
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Litvinova & Nadezhin (1985) caiibration Popov (1993) caiibration 




Figure 4. Explosion ener gy Eexp as a function of eject a mass Mgj. T he colors of th e individual supernovae are the same as in the right panel of Figure[T] We 
use the scaling relations of ILitvinova & NadezhinI <19851 left panel) and lPopo^ <19931 right panel). 




Figure 5. Plateau durations and expansion velocities from the theoretical models and the observations. The results of the theoretical non-rotating progenitor 
explosions of IDessart et al.l j2Q10l) are sho\vn with plus signs, stars, diamonds and triangles for explosion energies of 0.1, 0.3, 1.0 and 3.0 x 10^^ ergs. The color 
of the symbols indicate the ejecta mass. The obseiwations ai‘e shown with gray solid circles and their convex hull is marked with the gray polygon. The typical 
parameters of Type Il-Plateau supemovae (£p ^ 120 days, v ^ 3500 km s~^) can be explained by a range of models going from low Eexp and low Mej to 
moderate Eexp and high Mej. Finer details of the light curves and spectra need to be analyzed to break this degeneracy. 


energy (Eexp = 1.0 x 10®^ ergs) and relatively high ejecta 
mass (Mej > 10 Mq). In other words, very different sets of 
theoretical parameters will yield very similar observables, at 
least for tp and v (and likely also Lpi due to the observed 
Lpi-v correlation), implying that also the more sophisticated 
models are prone to the same degeneracy as the linear scaling 
relations. 


5. DISCUSSIONS AND CONCLUSIONS 

We show that the correlated uncertainties between param¬ 
eters derived from supernova light curves and velocities are 
strong and are oriented along the parameter correlations. The 
covariances between the quantities arise either due to uncer¬ 
tainties in the distance affecting two quantities in the same 
way (Lpi and Mni, Fig. [TJ, or due to degeneracies inher¬ 
ent to the physics of the supernova light curves (i?exp and 
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Mej, Fig. As a result, the statistical significance of these 
correlations is reduced, but the correlations cannot be fully 
explained away by the covariances. The correlation of Mej 
and Eexp is weaker than the other two correlations we investi¬ 
gated (Lp\, Eexp, and Mni) and its significance is sensitive to 
whether SN1980K (Type II-Linear) and SN1995ad (Type II- 
Plateau with short plateau duration) are included. Conversely, 
properly characterizing the uncertainty ellipsoids reveals an 
intrinsic width to the Lpi-MNi and Egy^p-Myn correlations 
(Figs. [U and |2]l. Now we explore the astrophysical implica¬ 
tions of these findings. 

In studies of the neutrino mechanism, most of Sexp comes 
from the neutrino-driven wind emanating from the nascent 
proto - neutron star (e.g. [Schecket'an 120061: lUgliano et alJ 
1201 21: iPeicha & ThompsonI 120151) . In a simple spherical 
picture, the evolution of the neutrino-driven wind is deter¬ 
mined by the thermo dynamic structure of the l ayers below 
the ejecta mass cut (IPeicha & Thomt)s'onll2015t) , while the 
ejected mass of ^®Ni depends primarily on the mass of the 
shock-heated ejecta exposed to sufficiently high temperatures 
(IWeaver & Woc)sl^ll980t IWooslevlll98^ iThielernann et alJ 
119901) . The intrinsic scatter in the E^xp-Myn relation there¬ 
fore implies that the progenitor structure below and above 
the mass cut cannot be fully des cribed by a single pa 


ramet e r, such as the compac tn ess (IO’Conn pr.&_Otl 


201 31: In akamura et al.l 120141: IPeicha & Th ompson 


2011 


2015 


Perego et al.ll2015l) . Using the results of IPeicha & ThompsCT 


dMl . we find a width of O.lOtg'oidex in their Ee^p- 
Mni correlation, which is consistent with the results pre¬ 
sented here. The theoretical pre dictions of the correlation 
slope (IPeicha & ThompsonI 120151) approximately agree with 
the observations pr esented here, b ut disagree with the infer¬ 
ences presented by iHamu'^ (120031) . A similar conclusion on 
the presence of multiple parameters determining the proper- 
ties of core-collapse supernova light curves was reached by 
ISanders et al.l (120151) . who quantified a dispersion in the de¬ 
cline rate-peak magnitude relation of hydrogen-rich super¬ 
novae. 

Changes in the intrinsic width of the E^xp-M^i correla¬ 
tion as a function of E^xp or other parameters can further 
constrain the explosion physics. For example, the spread in 
Mni at constant Eg^p could be caused by a varying amount 
of ®®Ni fallback in different progenitors. We would expect 
that the fallback will be generally less important at higher 
Eexp and the spread of Mni should thus increase as E^^p 
decreases. In principle, this is testable given a large set of 
well-observed supernova explosions. New un biased surveys 
of bright nearby supernova e such as ASAS-SN (iHoloien et al.l 
I20l4t1^appee et al.ll2014l) are particularly useful due to fea¬ 
sibility of detailed follow-up observations and the exploration 
of new parts of the parameter space (e.g., low-metallicity stel¬ 
lar environments). _ 

Contrary to the common picture (e .g. |Hegeretal .l |2003l: 
iNomoto et al.ll20d5 lUtrobin & Chugaill2009ll2014l) . there is 
little evidence from the parameterized studies of the neu¬ 
trino mechanism that the supernova properties such as Eg^p 
or Mni will strongly correlate with the mass of the progenitor 

(e.g.lO’ Connor & Ottl201 U lUgliano et al.l201^lBruenn et akl 


201 4: Nakamura et akl 2 01^ Pacha & ThompsonI 120151: 

Perego et al.l 120151: lErtl et al.l 120151) . because the ultimate 
fate of the star and the initiation of the explosion is set 
by the physics and the thermodynamic structure on the in- 
2.5 Mq of the progenitor, which is not monotonic 


ner 


with the initial mass met allicity or final hydrogen mass 
(ISukhbold & Wooslevll2^ . 

The degeneracy between Mej and Eg^p can be reduced by 
modeling finer features in the light curve and spectra, such as 
the O I 6303-6363 A line advocated bv iDessart et akl (120101) . 
However, characterizing these features generally requires sub¬ 
stantial observational effort on nearby objects (for example 
the 016303-6363 A line fully develops only ~ 300 days after 
explosion; iDessart et al.ll2010l) , along with confidence in the 
underlying physics and considerable effort in its numerical 
implementation. Surveys focusing on bright, nearby super- 
novae such as ASAS-SN (IHoloien et al.ll20l4lShat)pee et al.l 
12014h could be of help in this regard. It will be much harder 
to obtain such detailed information in the future dominated by 
primarily photometric discovery machines such as th e Large 
Synoptic Survey Telescope (e.g. Ilvezic et^ l2014ft . As a 
result, reliable recovery of supernova parameters from light 
curves and (potentially scarce) expansion velocities will only 
grow in importance, mandating rigorous uncertainty analysis 
that consistently include all relevant contributions to the final 
uncertainty budget, as we have done here. Detailed investi¬ 
gations of the supernova parameter covariances, degeneracies 
and multiple solutions based on detailed radiation hydrody¬ 
namic models is the next logical step in the preparation for 
the upcoming surveys. The understanding of the intricacies 
of the parameter recovery will not only yield greater physical 
understanding of the supernova population, but can influence 
back the design and strategy of the surveys currently in prepa¬ 
ration. 


ACKNOWLEDGEMENTS 

We thank Jujia Zhang for providing us with the velocities 
of SN2013am. We thank Chris Kochanek, Michael Strauss, 
Bruce Draine, and Todd Thompson for a detailed reading of 
the manuscript. We thank the referee for comments. OP ap¬ 
preciates the discussions with Adam Burrows. We thank the 
referee for comments that helped to improve the manuscript. 
Support for OP was provided by NASA through Hubble Fel¬ 
lowship grant HST-HF-51327.01-A awarded by the Space 
Telescope Science Institute, which is operated by the Asso¬ 
ciation of Universities for Research in Astronomy, Inc., for 
NASA, under contract NAS 5-26555. Support for JLP is 
provided in part by the Ministry of Economy, Development, 
and Tourism’s Millennium Science Initiative through grant 
IC120009, awarded to The Millennium Institute of Astro¬ 
physics, MAS. 


REFERENCES 


Anderson, J. R, Gonzalez-Gaitan, S., Hamuy, M., et al. 2014, ApJ, 786, 67 
Arcavi, I., Gal-Yam, A., Cenlco, S. B., et al. 2012, ApJ, 756, LL30 
Arnett, W. D. 1980, ApJ, 237, 541 

Bose, S., Kumar, B., Sutaria, F., et al. 2013, MNRAS, 433, 1871 
Bruenn, S. W., Lentz, E. J., Hix, W. R., et al. 2014, arXiv: 1409.5779 
Burrows, A. 2013, Reviews of Modem Physics, 85, 245 
Chatzopoulos, E., Wheeler, J. C., & Vinko, J. 2012, ApJ, 746, 121 
Dessart, L., Livne, E., & Waldman, R. 2010, MNRAS, 408, 827 
Dessart, L., Hillier, D. J., Waldman, R., & Livne, E. 2013, MNRAS, 433, 
1745 

Elmhamdi, A., Danziger, I. J., Chugai, N., et al. 2003, MNRAS, 338, 939 
Ertl, T, Janka, H.-T., Woosley, S. E., Sukhbold, T, & Ugliano, M. 2015, 
arXiv:1503.07522 

Faran, T., Poznanski, D., Filippenko, A. V., et al. 2014, MNRAS, 445, 554 
Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 
125, 306 

Hamuy, M. 2003, ApJ, 582, 905 





















































































8 


Pejcha & Prieto 


Heger, A., Fryer, C. L., Woosley, S. E., Langer, N., & Hartmann, D. H. 

2003, ApJ, 591, 288 

Hendry, M. A., Smartt, S. J., Crockett, R. M., et al. 2006, MNRAS, 369, 
1303 

Hogg, D. W., Bovy, J., & Lang, D. 2010, arXiv: 1008.4686 
Holoien, T. W.-S., Prieto, J. L., Pejcha, 0., et al. 2014, arXiv:1411.3322 
Insen-a, C., Pastorello, A., Turatto, M., et al. 2013, A&A, 555, AA142 
Ivezic, Z., Tyson, J. A., Abel, B., et al. 2014, arXiv:0805.2366 
leffreys, H. 1983, The Theoiy of Probability, 3rd ed., Oxford University 
Press 

Kasen, D., & Woosley, S. E. 2009, ApJ, 703, 2205 
Litvinova, 1.1., & Nadezhin, D. K. 1983, Ap&SS, 89, 89 
Litvinova, 1. Y., & Nadezhin, D. K. 1985, Soviet Astronomy Letters, 11, 145 
Nagy, A. P, Ordasi, A., Vinko, J., & Wheeler, J. C. 2014, A&A, 571, AA77 
Nakamura, K., Takiwaki, T, Kuroda, T, & Kotake, K. 2014, 
arXiv: 1406.2415 

Nomoto, K., Tominaga, N., Umeda, H., Kobayashi, C., & Maeda, K. 2006, 
Nuclear Physics A, 777, 424 
O’Connor, E., & Ott, C. D. 2011, ApJ, 730, 70 
O’Connor, E., & Ott, C. D. 2013, ApJ, 762, 126 
Pastorello, A., Zampieri, L., Turatto, M., et al. 2004, MNRAS, 347, 74 
Pejcha, O., & Kochanek, C. S. 2012, ApJ, 748, 107 
Pejcha, O., Thompson, T. A., & Kochanek, C. S. 2012a, MNRAS, 424, 1570 
Pejcha, O., Dasgupta, B., & Thompson, T. A. 2012b, MNRAS, 425, 1083 
Pejcha, 0., & Prieto, J. L. 2015, ApJ, 799, 215 
Pejcha, 0., & Thompson, T. A. 2012, ApJ, 746, 106 
Pejcha, 0., & Thompson, T. A. 2015, ApJ, 801, 90 


Perego, A., Hempel, M., Frohlich, C., et al. 2015, arXiv:1501.02845 
Popov, D. V. 1993, ApJ, 414, 712 
Poznanski, D. 2013, MNRAS, 436, 3224 

Press, W. H., Teukolsky, S. A., Vetterling, W. T, & Flannery, B. P. 1992, 
Cambridge: University Press, 2nd ed. 

Prieto, J. L., Stanek, K. Z., Kochanek, C. S., et al. 2008a, ApJ, 673, L59 
Prieto, J. L., Stanek, K. Z., & Beacom, J. F. 2008b, ApJ, 673, 999 
Prieto, J. L., Kistler, M. D., Thompson, T. A., et al. 2008c, ApJ, 681, L9 
Prieto, J. L., Lee, J. C., Drake, A. J., et al. 2012, ApJ, 745, 70 
Prieto, J. L., Brimacombe, J., Drake, A. J., & Howerton, S. 2013, ApJ, 763, 
L27 

Sanders, N. E., Soderberg, A. M., Gezari, S., et al. 2015, ApJ, 799, 208 
Scheck, L., Kifonidis, K., Janka, H.-T, Muller, E. 2006, A&A, 457, 963 
Shappee, B. J., Prieto, J. L., Grupe, D., et al. 2014, ApJ, 788, 48 
Smartt, S. J. 2009, ARA&A, 47, 63 

Spiro, S., Pastorello, A., Pumo, M. L., et al. 2014, MNRAS, 439, 2873 
Sukhbold, T, & Woosley, S. E. 2014, ApJ, 783, 10 
Thielemann, F.-K., Hashimoto, M.-A., & Nomoto, K. 1990, ApJ, 349, 222 
Ugliano, M., Janka, H.-T, Marek, A., & Arcones, A. 2012, ApJ, 757, 69 
Utrobin, V. P, & Chugai, N. N. 2009, A&A, 506, 829 
Utrobin, V. P, & Chugai, N. N. 2014, arXiv:1411.6480 
Tremaine, S., Gebhardt, K., Bender, R., et al. 2002, ApJ, 574, 740 
Weaver, T. A., & Woosley, S. E. 1980, Ninth Texas Symposium on 
Relativistic Astrophysics, 336, 335 
Woosley, S. E. 1988, ApJ, 330, 218 

Zampieri, L., Pastorello, A., Turatto, M., et al. 2003, MNRAS, 338, 711 
Zhang, J., Wang, X., Mazzali, P. A., et al. 2014, ApJ, 797, 5 


